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f-r—\  Introduction 

<'in  a  recent  series  of  papers  ([2],  [3],  [4])  the  author 
has  investigated  the  low  frequency  inverse  scattering  problem 
for  "hard"  and  "soft"  infinite  cylinders  and  "soft”  obstacles 
in  space.  It  is  the  purpose  of  this  paper  to  complete  the 
story  and  consider  the  inverse  scattering  problem  in  space 
where  the  obstacle  is  "hard",  i.e.  the  boundary  condition  is 
of  Neumann  type.  A  basic  result  of  the  analysis  to  date  is 
that  the  low  frequency  inverse  scattering  problem  can  be 
stabilized.  In  the  case  of  domains  in  the  plane  this  is 
accomplished  by  transmitting  two  plane  waves  from  different 
directions,  where  the  second  plane  wave  is  used  to  determine 
the  transfinite  diameter  of  the  obstacle# (  [2] ,  [3]).  Due  to 
the  lack  of  the  availability  of  conforrfial^  mapping  techniques, 
this  approach  does  not  extend  to  the  three  dimensional  inverse 
scattering  problem  and  alternate  methods  must  be  developed. 

In  the  case  of  a  "soft"  obstacle  (i.e.  Dirichlet  boundary  data) 
stability  was  obtained  by  assuming  an  "a  priori"  knowledge  of 
the  radius  of  a  ball  containing  the  unknown  scattering  obstacle 
in  its  interior  and  using  the  fact  that  the  boundary  of  the 
scattering  obstacle  is  a  level  surface  of  the  limiting  static 
potential  problem  ((4]).  Since  the  static  solution  is  simply 
the  conductor  potential  of  the  obstacle  irregardless  of  which 
direction  the  incoming  wave  impinges  upon  the  obstacle,  no  new 
information  is  gathered  by  sending  in  plane  waves  from  different 
directions  (however  see  Section  IV  of  this  paper) .  In  the 
present  paper  we  shall  be  considering  the  inverse  scattering 
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problem  for  a  "hard"  obstacle  in  space  and  as  will  be  seen  the 
situation  is  considerably  altered.  In  particular,  in  order  to 
achieve  stability  it  is  necessary  to  not  only  assume  an  "a 
priori"  knowledge  of  the  radius  of  a  ball  containing  the 
scattering  obstacle  in  its  interior,  but  also  on  "a  priori" 
bound  on  the  gradient  of  the  velocity  potential  of  the  total 
field  for  small  values  of  the  wave  number.  Further  complications 
arise  from  the  facts  that  the  boundary  of  the  obstacle  is  no 
longer  a  level  surface  of  an  appropirate  potential  problem  and 
the  solution  of  the  limiting  static  potential  problem  is  simply 
u_(x)=l  which  gives  no  information  on  the  shaoe  of  the  obstacle. 
These  problems  will  be  circumvented  by  identifying  the  second 
term  in  the  low  frequency  expansion  of  the  total  field  as  the 
velocity  potential  of  an  incompressible,  irrotational  fluid 
flow  past  the  scattering  obstacle  and  using  this  ^function  to 
determine  the  streamlines  of  the  flow  outside  the  given  ball 
containing  the  obstacle  in  its  interior.  In  this  case  sending 
in  plane  waves  from  different  directions  is  equivalent  to  fluid 
flowing  past  the  obstacle  in  different  directions,  and  knowledge 
of  the  resulting  streamlines  then  gives  more  information  on  the 
shape  of  the  obstacle.  In  order  to  estimate  the  actual  shape 
of  the  obstacle  an  important  role  is  played  by  the  fact  that 
from  the  far  field  data  it  is  possible  to  deduce  the  volume  of 
the  obstacle.  Hence  from  a  knowledge  of  the  streamlines  and 
the  volume  of  the  obstacle,  accurate  estimates  can  be  made  on 
its  shape. 
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i  II .  Reduction  of  the  Inverse  Scattering  Problem 

;  to  a  Nonlinear  Moment  Problem 

The  aim  of  this  section  is  to  reduce  the  inverse  scattering 

‘  i 

problem  for  a  "hard"  obstacle  to  a  nonlinear  moment  problem 

1  involving  the  velocity  potential  of  an  incompressible,  irrota- 

| 

tional  fluid  flow  past  the  unknown  scattering  obstacle.  We 
shall  assume  that  the  scattering  obstacle  D  is  a  bounded  simply 

f 

'  connected  domain  containing  the  origin  with  smooth  boundary  SD. 

In  order  to  achieve  our  aim  it  is  necessary  to  first  consider 
the  direct  scattering  problem,  i.e.  to  find  a  solution 
u  £  C2  (R3\D)  n  C1  (R3\D)  such  that 


is  well  known  (c.f.  [9])  and  our  first  aim  is  to  find  a  low 
frequency  approximation  to  u  evaluated  on  3D.  To  this  end  we 


note  that  from  [1]  and  [7]  we  can  reformulate  (2.1a)  -  (2. Id) 
as  the  integral  equation. 


ui(x) 


♦  f i“ 


(y) 

*  v 

3D 


u(x)  ]  4”  — - — r  ds(y) 
^  jx-yj 

X  X 


~  f  U(v)  -I— 
4  IT  J  ^  <3\> 

3D 


[exp  [ik  !  £-y !  ]  -l] 


(2.2) 


ds(y)  =  u(x) 


x  -  y| 


for  x  £  R2\D ,  or  in  obvious  operator  notation 

% 

ux  +  L,  [u]  =  u.  (2.3) 

\ 

The  important  point  here  is  that  the  spectral  radius  of  L  is 
less  than  one  and  that 

M  -  L0  | |  =  0(k2)  (2.4) 

where  j  !  •  ! |  denotes  the  maximum  operator  norm  ([7]) .  Hence 
for  k  sufficiently  small,  (2.3)  can  be  solved  by  successive 
approximations  and  I!  (I  -  £k) _1  -  (I  -  Lg)"1  i!  -  0(k2)([8], 
p.  164).  In  particular  for  u1(x)  »  elkxl  x  *  (x.  ,  x,,  x,) , 

v  ±  4  j 

and  k  sufficiently  small,  we  have  for  xe3D  that 


-5- 


u(x)  > 

x  n*0  ^ 


2  tl  +  i  kxxl  + 


0(k2) 


[1  +  i  kxx]  +  0(k‘) 


(2.5) 


n=0 


1  +  ik  u  (x)  +  0  (k  ) 
°  ^ 


where  can  be  identified  as  the  uniaue  harmonic  function 
o 

defined  in  the  exterior  of  D  such  that 


u  (x) 
°  <x. 


*1  +  uo 


(x) 


in 


r3\d 


A3  u0  «  0 

M  <*>  -  0 

3 

lim  uQ  (x) 
r-<=° 


in  R3\D 


for  xc  3D 

0, 


(2.6a) 

(2.6b) 

(2.6c) 

( 2 . 6d) 


i.e.  uq  is  the  velocity  potential  for  an  incompressible,  irrota- 
tional  fluid  flow  past  D  having  constant  velocity  in  the  x^ 
direction  at  infinity. 

We  now  proceed  to  the  reduction  of  the  inverse  scattering 
problem  to  a  nonlinear  moment  problem.  From  (2.2)  we  have  that 
((51,  ?.  316) 


u 


s 


(X) 

V 


~ex?  lifc  ^  (y) 

'  x  -  y ;  J 


T3T(gW 


; 


F(n;k)  +  0  (i2) ?  n 

-v  'v, 


(2.7) 


where 


F  ( n ;  k )  =  j-  f  u (y)  1—  exp  [-ikn'y]  ds(y) 


(2.8) 


The  inverse  scattering  problem  we  are  considering  is  to  determine 
D  from  a  knowledge  of  the  far  field  pattern  F  for  all  directions 
n  and  low  values  of  the  wave  number  k.  To  this  end  we  write  n 
and  y  in  spherical  coordinates  as 


n  =  (sin  9  cos  sin  3  sin  o,  cos 3) 

^  (2.9) 

y  =  p(sin  9'  cos  $  ’ ,  sin  9'  sin  cos  9') 


and  expand  F  in  a  series  of  spherical  harmonics 


'*>  -  2  2  s 


,,  v  _m  ,  , .  im® 

„„(k)  p„  (cos  9)  e 
mn  n 


n*0  m=-n 


(2.10) 


Then  we  can  easily  deduce  from  (2.8),  (2.10)  and  the  relation 
(c.f.  [6]) 
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e'il!n  */jf?  2  <-i>n  ,kpt  pn(MsY) 

n=o  2 

(2.11) 

n  •  y  =  p  COS  Y 

v  a, 

that 


a 

nm 


(k) 


(kp)  • 


_m,  .  , ,  -imp ' 

P  (cos9')e 
n 


ds  (y) 


(2.12) 


where  we  have  used  the  standard  notation  for  associated  Legendre 
polynomials,  Legendre  polynomials,  and  Bessel  functions.  Since 


—  i  (kp)  *  — i -  +  0  (kp)nT,i  (2.13) 

✓Jc?  n+  T  2n+T  r(n+A) 


we  have 


anm(k) 

nm 


from  (2.12)  and  (2.13)  that  for 

i~n/^ntn 

2n+1  F (n+T) 

+  0  (kn+2) 


3D 


v 


m 


n>l 


(cos9 ’ 


ds  (y) 

(2.14) 


and  since  the  quantity  in  brackets  is  harmonic,  (2.5)  and  (2.14) 
imply  that 


mm 
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a  (k)  = 
nm 


i-n+1/¥~kn+l 
2n+1  P(n+f)  3D 


I  u„(y 


+  0(k 


o  r  av 


n+2 ) 


pn  P®  (cosS ' ) e~im^ '  • 


or  n>_l.  In  the  case  n=0  a  similiar  analysis  shows  that 

/  1_ 

J  3v 


a  (k)  " 
oo 


/t  k" 


ar  (4-)  3d 

If  we  now  define 


o2  ds  +  0  (k3) 


nm 


oo 


lim 

k-o 


2n+1r (n+  V) 

i*n+1/7kn+1  “nm 


a _ (k) 


;  n>  1 


lim 

k-*-o 


4r(f) 

3k2 


a  (k) 
oo 


nm 


we  can  rewrite 

(2 

.15) 

in  the  form 

-  J*  «0<y> 

3__ 

3  v 

n 

P 

Dm 

*  n 

(cos9  ' )  e  1 

3D  % 

m 

J 

ds  (y)  ;  n>l 


and  note  that  (2.16)  yields 


oo 


-*/ 

3D 

.iff- 

D 

■  //■ 


3  _  2 


p  ds 


l3P2dv 


ds  (y) 
(2.15) 


(2.16) 


(2.17) 


(2.18) 


(2.19) 


vo 1 ume  of  D . 


Hence  we  have  reduced  our  inverse  scattering  problem  to 


that  of  determining  D  from  the  nonlinear  moment  problem  (2.18) 
along  with  the  fact  that  we  now  know  the  volume  of  D  from  the 
far  field  data  y  .  Note  that  although  the  ynm  are  explicitly 
computable  from  the  far  field  pattern  F,  small  errors  in 
measuring  the  coefficients  a  (k)  will  result  in  large  errors 
in  the  numbers  ymn  if  n  is  large.  Hence  in  practice  we  must 
assume  that  only  a  finite  number  of  the  unm  are  known  and  for 
the  sake  of  simplicity  we  shall  assume  that  these  are  known 
exactly.  Observe  that  the  "moment"  problem  (2.13)  is  nonlinear 
since  both  uQ  and  the  region  of  integration  depend  on  D. 

Ill .  Analysis  of  the  Nonlinear  Moment  Problem 

We  shall  now  show  how  information  on  the  shape  of  D  can 

be  extracted  from  the  nonlinear  moment  problem  (2.18).  First 

supoose  that  the  moments  y  are  known  for  n  *  1,  2...N, 

nm 

-n<m<_n,  and  that  constants  a  and  M  are  known  such  that 

(1)  D  =  B  =  {x  :  i  x!<a} 

■v  'V, 

(2)  max  !  grad  u  !  <kM  for  0<k<k 

3  ~  - ° 

R  \  3 

where  k  is  a  fixed  constant  and  M  is  independent  of  k.  Note 
o 

that  from  (2.5)  we  have  from  the  second  condition  that 


max 


•  grad  u  i  <  M. 


(3.1 
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For  xeSB  we  have  that 
U  (x)  =  K1  + 


u  has  an  expansion  of  the  form 
o 

oo  n 


2  2 

n=0  m=-n 


,  ini 

b_o  P  (cos0)e 
run  n 


(3.2) 


Furthermore  from  Green's  formula,  (2.6c)  and  (2.18)  we  have  that 
for  n^l 


nm 


{u  ~  [pVtcossMe'1^'  ] 
O  a  o  n 

3u 


(3.3) 


^-[pnP:n  (cose  1  )  e  ]}ds 
ap  n 


and  from  (3.2),  (3.3)  and  the  orthogonality  of  the  associated 
Legendre  polynomials  we  can  conclude  that 


unm  =  4"  bnm  ?  * 


(3.4) 


For  n=0  we  have  from  Green's  formula  that 

I*  3u„ 


0  =  - 


3p 


3B 


ds 


0  0 


=  b 


oo 


*  4j  b 


inededi 


oo  ’ 


(3.5) 


nm 


i.e.,  bQO  =  0.  Hence  since  the  a _  are  assumed  known  for 

n=l, 2, . . . ,N, -n<m<n,  we  know 

N  n 


U^(x)  =  X1  +  2  2 


n=0  m=-n 


,  -n-l_m,  , .  imo 

bnm3  Pa(cos8)e 


(3.6) 
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and  from  assumptions  1)  and  2)  above  we  can  easily  obtain  explicit 

N 

L-,  error  estimates  for  grad  (uQ-u0)  in  terms  of  the  constants  a, 

3  N 

M  and  N.  In  particular  if  x=(x.  ,x,,x,)  then  t— (u_-u‘ )  ,  i*l,2,3, 

x  3  ®  o 

is  harmonic  and  has  an  expansion  of  the  form 


2  2 

n=N+l  m=-n 


,(i)a-n-lpm  imo 

nm  n 


(3.7) 


Hence  we  can  conclude  by  a  short  calcuation  (c.f.  [4])  that 


4rp  |  2C |  — p 


J  | grad (u„-u^) j 2ds 


o  o 


OD 

<  6oM  X 


n=N+l 


( 2n+l) (— ) 2n+2  (3.8) 

p 


and  thus  we  can  decermine  how  many  Fourier  coefficients  of  the  far 
field  pattern  are  needed  in  order  to  approximate  grad  uQ  for 
|x|>a  >a  to  obtain  a  given  degree  of  accuracy  in  the  L2  sense. 

Identifying  uQ  as  the  velocity  potential  of  an  incompressible, 
irrotational  fluid  flow  past  D,  we  now  note  that  the  unknown 
obstacle  D  is  a  surface  containing  streamlines  of  the  fluid  flow 
past  D.  Hence  if  we  can  determine  streamlines  lying  in  the  exterior 
of  D  this  will  give  us  information  ou  the  shape  of  D.  To  this  end 
let  xQ  be  a  point  in  S  \B.  Then  the  streamline  passing  through  xQ 
can  be  determined  by  solving  the  system  of  ordinary  differential 
equations 

dS  . 

dt  =  grad  U 

x (0 )  =  xQ 


(3.9) 


-12- 

where  x=x(t)  is  a  parametric  representation  of  the  given  stream¬ 
line.  Since  grad  u  is  only  known  to  a  given  degree  of  accuracy 
for  x  in  £3\B,  (3.9)  will  only  allow  us  to  construct  the  stream¬ 
lines  lying  in  the  exterior  of  B.  However,  we  do  know  the 
volume  of  D  from  (2.19),  and  hence  by  extrapolation  we  can  deter¬ 
mine  an  approximation  to  the  shape  of  D.  Note  that  by  sending  in 
plane  waves  from  different  directions  we  can  also  determine  the 
streamlines  corresponding  to  a  fluid  flowing  past  D  from  different 
directions  and  this  can  be  used  to  refine  the  initial  estimate  for 
the  shape  of  D. 

IV.  Passing  Remarks  on  the  Dirichlet  Problem 

We  noted  in  the  Introduction  that  for  the  Dirichlet  problem 
the  solution  of  the  static  potential  problem  corresponding  to  a 
single  incoming  plane  wave  is  the  conductor  potential  and  hence, 
in  contrast  to  the  Neumann  problem,  no  new  information  is  gained 
by  sending  in  plane  waves  from  different  directions.  This  parti¬ 
cular  problem  can  be  overcome  if  two  plane  waves  are  transmitted 
from  opposite  directions,  e.g.  e1**!  and  e~^'<xl,  and  the  corre¬ 
sponding  solutions  of  the  scattering  problems  added  together  in 
such  a  way  that  the  combined  incident  field  is  i-sin  kxj_.  The 
analysis  of  [4]  can  then  be  used  to  show  that  the  corresponding 
static  potential  problem  or.  this  case  is 


S  3 

u  (x)  -  x,  +  u  (x)  in  1  \D 

O  -  i  O' 

(4.1a) 

43uo  *  0  in  S3  D 

(4.1b) 

u  (x)  =  0  for  x£jD 

(4.1c) 

lini  u3(;c)  =  0 

-*-oo 

(4. Id) 

The  boundary  of  the  obstacle  is  again  a  level  surface  and  except 
for  this  boundary  the  level  surfaces  change  if  the  plane  waves 
are  sent  in  from  different  directions  than  along  the  axis. 
Hence  if  this  is  done  more  information  can  be  obtained  on  the 
shape  of  D. 
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